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ABSTRACT 

The magnetic flare avalanches model of Poutanen & Fabian for X-ray variability 
of accreting black holes is combined with computations of vertical structure of illu- 
minated accretion discs in hydrostatic equilibrium. The latter predict the existence of 
a hot ionized skin, due to the thermal instability of X-ray illuminated plasma. The 
presence of such ionized skin, with properties dependent on disc radius, introduces 
a dependence of the emitted X-ray spectrum on the position on the disc. If the po- 
sition is related to the time scale of the flares, the X-ray energy spectra (both the 
primary continuum and the reprocessed component) gain an additional dependence 
on Fourier frequency, beside that resulting from spectral evolution during a flare. We 
compute the Fourier frequency resolved spectra in this model and demonstrate that 
the presence of the hot skin introduces trends opposite to those observed in black hole 
binaries. Furthermore, the flare profile is strongly constrained, if the Fourier frequency 
dependence due to spectral evolution is to agree with observations. 
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1 INTRODUCTION 

X-ray emission from accreting compact objects carries in- 
formation about geometry and physical conditions in the 
immediate vicinity of central objects. Unfortunately, the in- 
formation is "encrypted" into a series of numbers, each num- 
ber (event) corresponding to a detection of a photon from 
the source. Two main characteristics of each event are the 
photon energy and arrival time (we neglect here polariza- 
tion of radiation). One way of analyzing the data stream is 
then to project the events onto the energy axis, thus sup- 
pressing the timing information. Thus obtained energy spec- 
tra give a first, rough insight into physical mechanism of 
emission of radiation. For example, observations of cutoffs 
at ~ 100 keV in low/hard state of accreting black holes 
demonstrate that the Comptonizing electrons have a ther- 
mal rather than non-thermal distribution. From modelling 
the spectra a mean plasma temperature and scattering op- 
tical depth can be inferred. However, any possible distribu- 
tions of plasma temperature and density cannot be recovered 
from time-averaged spectra (e.g. Kazanas, Hua & Titarchuk 
1997; Coppi 1999). 
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In an analogous way, projecting the data onto the time 
axis produces a light curve, which can be analyzed either 
in the time or Fourier domain. The results, too, enable de- 
termination of certain average properties. For example from 
power spectral density (PSD), the root mean square vari- 
ability or constraints on distribution of shots' time scales in 
the shot noise model can be found (Lehto 1989; Lochner, 
Swank & Szymkowiak 1991). 

Further progress can be made by including at least lim- 
ited information from along the other parameter axis. For 
example, cross-correlation analysis (both in time and Fourier 
domain) reveal that higher energy X-rays lag the lower en- 
ergy ones (Miyamoto & Kitamoto 1989), which one can at- 
tempt to interpret as e.g. spectral evolution during a shot 
(flare). Then the energy dependent auto-correlation func- 
tions give average shot profiles as a function of energy (see 
review in Poutanen 2001 and references therein). 

Similarly, studying a sequence of energy spectra, each 
one averaged over a certain short time interval, enables de- 
termination of not only the mean spectral properties, but 
also "higher orders" statistics properties. The most impor- 
tant result that has emerged recently is that spectral prop- 
erties of a given source vary with time, but correlations be- 
tween parameters exist. The power law slope, F, is corre- 
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lated with the amphtude of reflection R and the width of 
the Fe Kc hne (Zdziarski, Lubihski & Smith 1999; Gilfanov, 
Churazov & Revnivtsev 1999; Revnivtsev, Gilfanov & Chu- 
razov 2001; Zycki, Done & Smith 1998). Interestingly, the 
same correlations seem to exist between sources in a sample 
(Zdziarski et al. 1999; Lubiiiski & Zdziarski 2001). 

A number of geometrical/physical scenarios can be built 
based on the R-T correlation alone. The crucial element of 
all models is the feedback between the hot. X-ray emitting 
plasma and the cooler, optically thick plasma: most of the 
soft photons crossing the hot plasma have to come from 
thermalization/reprocessing of the illuminating hard X-rays 
(Poutanen 1998 and references therein). Given this condition 
alone, models can be constructed belonging to two main 
classes, (1) with the standard accretion disc disrupted at 
a certain radius, and (2) with the standard accretion disc 
extending to the last stable orbit: 

• standard accretion disc disrupted at a given radius, 
and filled with hot, optically thin, X-ray producing flow 
(Esin, McCUntock & Narayan 1997; Poutanen, KroUk & 
Ryde 1997; Rozariska & Czerny 2000), with the amount of 
overlap between the two phases controlling the correlations, 

• hot, optically thin, quasi-spherical flow with cold 
cloudlets inside (Rees 1987; Kuncic, Celotti & Rees 1997), 
with the cloudlets flUing factor controlling the correlations, 

• spherical accretion of cold clouds onto a central X-ray 
source ( Collin- Souffrin et al. 1996; Malzac 2001), with the 
covering factor of the clouds as the control parameter, 

• mildly relativistic outflow of the X-ray emitting plasma, 
above a standard accretion disc, with the outflow veloc- 
ity controlling the correlations (Beloborodov 1999a,b; Ja- 
niuk, Czerny & Zycki 2000; Malzac, Beloborodov, Poutanen 
2001), 

• hot, ionized skin on top of the accreting disc, due to 
thermal instability of the X-ray illuminated plasma, with 
the skin thickness as the control parameter (Rozahska & 
Czerny 1996; Nayakshin, Kazanas & Kallman 2000; Zycki 
& Rozahska 2001). 

Also, models for time variability predicting correct 
PSD, time lags and coherence function can be constructed 
for the two main geometries (Poutanen & Fabian 1999, 
hereafter PF99; Bottcher & Liang 1999; review in Pouta- 
nen 2001). Therefore, based solely on the limited spec- 
tral/temporal data distinguishing between these scenarios 
appears very difficult, if at all possible (see e.g. discussion 
and application to Cyg X-1 in Di Salvo et al. 2001). Employ- 
ing more sophisticated tools, involving both spectral and 
temporal properties is clearly required. 

The Fourier frequency resolved (/-resolved hereafter) 
energy spectra seem to be a potentially powerful tool to 
apply to the problem of distinguishing the models. The /- 
resolved spectra are energy spectra in the limited range of 
Fourier frequency and so, roughly speaking, represent spec- 
tra at different time scales of variability (Revnivtsev, Gil- 
fanov & Churazov 1999; Gilfanov, Churazov & Revnivtsev 
2000a,b). The /-resolved spectra found from data of X-ray 
binaries show a number of features and correlations. Most 
importantly, the .R-F correlation is fulfilled by the /-resolved 
spectra in the sense that the higher the Fourier frequency, 
the harder the spectrum, and the smaller the amplitude of 
reflection (Revnivtsev et al. 1999). This result indicates that 



the same mechanism which is responsible for the correlation 
on long timescales (~ 10^ sec ~ 10^-10^ tdyn, tdyu = 
being the dynamical timescale), acts also on much shorter 
time scales, at least down do ~ 0.01 sec, the time bin of the 
Revnivtsev et al. analysis. This fact provides the motivation 
for our attempt to connect the spectral and timing models 
and predict their combined results. Therefore in this paper 
we consider a simple combination of the model of PF99, 
with the structure of illuminated accretion discs in hydro- 
static equilibrium. The former is the most complete X-ray 
emission model so far, since it addresses both the spectral 
and timing properties of Comptonization emission. The lat- 
ter naturally predicts the appearance of a hot, ionized skin 
on top of the disc (Rozariska & Czerny 1996; Nayakshin et 
al. 2000), as a result of the thermal instability of the X- 
ray irradiated plasma (Field 1965; Krolik, McKee & Tarter 
1981). The ionized skin reduces the effectiveness of repro- 
cessing and thermalization of the hard X-rays. The flux of 
soft photons is reduced, which may lead to a harder con- 
tinuum spectrum. If the ionized skin is thickest close to 
the black hole (Nayakshin 2000; Zycki & Rozariska 2001), 
then the hardest spectra, with lowest reflection amplitude 
are expected to originate there. If the short time scale (high 
Fourier frequency, /) flares originate close to the black hole, 
then one may expect the high-/ spectra to be harder, and 
show less reflection compared to low-/ spectra, precisely as 
observed from accreting black hole binaries. 

The plan of the paper is as follows: in Section ^ we de- 
fine the model, in Section |^ we perform some test computa- 
tions to gain deeper insight into the properties of frequency- 
resolved spectra. Section ^ contains results of full model cal- 
culations and their comparison with data; the results are 
discussed in Section pi 



2 THE MODEL 
2.1 Time variability 

We envision that the X-rays are produced in flares occurring 
above an accretion disc. The variability part of our model is 
therefore identical to that of PF99 (see also Stern & Svens- 
son 1996), so only a brief description is given here for com- 
pleteness. There is a number of spontaneous flares, A per 
second, each of which can give rise to an avalanche of stim- 
ulated flares. The number of stimulated flares is drawn from 
a Poissonian distribution with the mean /x. Each stimulated 
flare is delayed from the stimulating flare by an interval At, 
drawn from a Poissonian distribution of mean ari , where ri 
is the time scale of the stimulating flare while a is a param- 
eter. The time scales of all flares are drawn according to a 
power law probability distribution cc t^"^ between Tmin and 

As in PF99 all flares have the same profile, which 
describes the heating of electrons as a function of 
time, Ihit)- Here, Zh is the compactness parameter, I = 
LaT / {Raa.reme<?), where L is the luminosity and -Rflarc is the 
radius of a fiare. The flare proflle turns out to be decisive 
for the obtained dependence of energy spectra on Fourier 
frequency, therefore we will employ two of prescriptions: 

• the function used by PF99, 

Zh,i(t) oc (t/r)^exp(-i/r), (1) 
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• a double exponential function defined as 



exp{t/Ti), for t < 
exp(— t/rd), for t > 0, 



(2) 



where Tr and represent the rise and decline time scales 
and Tr ^ Td (Maccarone, Coppi & Poutanen 2000). 

The functions are normalized to give an assumed value of 
peak compactness, /h,max, for which we adopt Zh.max = 10. 
The flare decline time scale is assumed Td — O.Iti. 

The flux of soft photons through the flare is given as 
ls{t) = ls,o + D(t)lh{t) (PF99), where the second term de- 
scribes the feedback due to reprocessing of the hard X-rays. 
Here, our first extension of the PF99 model is introduced; 
the feedback expression contains a factor describing the ef- 
fectiveness of thermalization due to the ionized skin, 



m = Sit) 



Do 



l + 3[H{t)Y 



(3) 



Here S{t) — ^[ThotC^)] is the fraction of the illuminating hard 
X-rays that get through the ionized skin to the cold disc, 
as a function of Thot (equation |l^), Do is a parameter, while 
H{t) is the height of the flare above the disc. 



2.2 Radial dependencies 

The crucial extension of the PF99 model which needs to be 
made here is to relate the time scale of a flare to its radial 
location on the accretion disc. We will consider two cases: 

(i) There is a unique correspondence between the time 
scale of an individual flare (both spontaneous and stimu- 
lated) and its location on the disc. Since we do not assume 
any correlation between time scales of a stimulating and 
stimulated flare, this means that spatial locations of fiares 
within an avalanche are not correlated. 

(ii) The unique correspondence between r and r holds 
for spontaneous flares, while all fiares stimulated by a given 
spontaneous flare occur at the location of the latter (Merloni 
& Fabian 2001). 

The essential assumption here is that there is at least partial 
correspondence between flare's time scale and its radial po- 
sition. To quantify the discussion we assume that the flare 
time scale has the same radial dependence as the inverse of 
the Keplerian frequency, r oc r^^^. Moreover, the shortest 
flare time scale is assumed to correspond to the marginally 
stable orbit at 6rg — (SGM/c? . This gives 



( — ) 



2/3 



(4) 



One consequence of such assumption is the existence of 
the maximum radius of the flares-covered disc, rmax = 

6rg (Tmax/Tmin)^''^ ~ 170rg. 

Some other parameters of the model need to be consid- 
ered as functions of radius, too. The compactness 

0"T 



TT-Rflarc-fdisc- 



(5) 



(see equation 1 in PF99) is a function of radius, since both 
the internal disc emission, Fdisc, and the flare's radius, i?flarc 
may in general be functions of radius. For the internal 
disc emission we adopt the standard, Newtonian formula 
(Shakura & Sunyaev 1973) 



-Fdisc (r) = e-Ftot(r) = C; 



GMM 




(6) 



where ^ is the fraction of total energy dissipated inside the 
disc (1 — C dissipated in the active corona). We will assume 
that the radius of a flare, -Rflaro, is constant with the radial 
position. The value of -Rflarc will be derived from global en- 
ergy budget for the hot plasma: the total, time-integrated 
luminosity of all flares equals energy dissipation in the hot 
phase. 



(7) 



where the sum is over all fiares, Ltot ~ f,°^ 2-Kr Ftotir) dr 

J org ^ ^ 

is the total accretion luminosity, while T is the total time 
interval. Note that the radial dependence of Zs,o introduces 
slight overall radial dependence of spectra, since the time 
evolution of electron temperature depends on the value of 



2.3 The energy spectrum 

2.3.1 Comptonized continuum 

We assume that the continuum spectrum is due to Comp- 
tonization of soft photons from the disc and most of these 
photons come from reprocessing of the hard X-ray flux. The 
Comptonized spectrum is computed using the thComp code 
of Zdziarski, Johnson & Magdziarz (1996), based on solution 
of the Kompaneets equation. This is parameterized by the 
photon index of the spectrum, F and plasma temperature, 
fcTc. These parameters are computed using formulae from 
Beloborodov (1999b): F = 2.33(^h/^s)"'''^ F = (9/4)y-2/^ 
and y = 4(6^0 -I- 46'g)Tcs(Tos + 1), where y is the Comptoniza- 
tion parameter and 6^ = kTc/nicC^. We assume Tos ~ 1.8 
(PF99). 



2.3.2 The reprocessed component 

From equation (1) of PF99 one can derive the flux of radi- 
ation illuminating the disc in the immediate vicinity of the 
flare. 



-Firr(t) 



/h(t) vDo 



UTTV ijflarc 1 + 3 [H {t)Y 



(8) 



where the height of the flare above of the disc is H(t) = 
Hot/r, with rj — 0.5 and Do and -ffo as parameters. From 
the known Firr the thickness of the ionized skin can be com- 
puted at the considered radius (0.25-Flrr is used in the actual 
computations, representing spatially average flux). This is in 
general a complicated task, requiring solving simultaneously 
equations of the vertical structure of the disc and the ionized 
skin (Rozariska & Czerny 1996; Rozariska 1999; Nayakshin 
et al. 2000; Zycki & Rozariska 2001). The computations can- 
not be performed in real time because of the large number of 
spectra needed. Simulating a 256 sec lightcurve with 1/128 
sec time resolution requires computing more than 10^ spec- 
tra. Therefore, we use the code of Zycki & Rozariska (2001) 
to generate a grid of Thot as a function of Firr and r. A fitting 
formula is then used to describe the results. 
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Figure 1. Thomson thickness of the hot slcin as a function of 
the illuminating flux and disk radius computed using the code of 
Zycki & Rozahska (2001). Fitting formula (Eq. ^) has been used 
as a representation of the results (curves) for a given value of the 
spectral index of illuminating radiation. Power law with F = 1.65 
was used in the presented case. 
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(9) 



erg cm ^ s ^ 



W = 1.05, and iV = 
The accuracy of the formula is 



where Fo = 10 
1.034 + 0.0937 log(r/2r, 
demonstrated in Fig. ^. 

Knowing rhot we estimate the effective amplitude of the 
cold reprocessed component as 



R 



[Ptr('rhot)]^ 



(10) 

1 +Prf(-rhot) 

where ptr(T"hot) is the probability of an illuminating photon 
being transmitted through the hot skin to the cold disc, 
while Prf(''"hot) = 1 — Ptr(''"hot) IS the probability of the pho- 
ton being back-scattered in the hot skin and contributing 
to the primary continuum (Zycki & Rozahska 2001). The 
square of pti in the numerator comes from the fact that the 
reflected photons have to go twice through the hot skin. The 
probabilities are estimated using a Monte Carlo code for a 
purely scattering atmosphere. The above prescription for R 
is obviously rather approximate, so we do not attempt to e.g. 
reproduce the correct observed amplitude of the reprocessed 
component. The function ^(rhot) describing the dependence 
of the soft flux from thermalization on rhot (see Eq. ^) is 
taken from Poutanen (2002), 

7-hexp(-10rhot) 



5'(rhot) = (l-a); 



(11) 



2[4-h3rhot(l-a)]' 
where a is the cold reflection albedo, a ~ 0.2. 

The spectrum of the reprocessed component is con- 
structed according to the simple formula of Lightman & 
White (1988), 



•S'rcfl(F) = - — ' 5prim(-E'), 6 



^abs ~\~ l^c; 



(12) 



where Kabs(£') and Kos(F) are the photo- absorption and elec- 
tron scattering opacities, respectively. "Cold" matter ab- 
sorption opacities and elements abundances from Morrison 
& McCammon (1983) are used. Eq. gives the angle aver- 
aged reflected spectrum and it is valid up to ~ 15 keV, i.e. 
as long as the electron scattering can be treated as elastic. 
Formula ( p^ ) is multiplied by a simple exponential cutoff 
to mimic the proper shape of the reflected continuum at 
higher energies. This rough approximation to the reflected 
spectrum is sufficient for our purposes, since the spectra only 
up to ~ 30 keV will be presented, and it is not feasible to 
use e.g. the proper angle-dependent functions of Magdziarz 
& Zdziarski (1995). The Fe Kc fine is added to the repro- 
cessed continuum. Its equivalent width (EW) is computed as 
a function of the spectral slope, F, based on computations by 
Zycki & Czerny (1994). We note a number of simplifications 
involved in the above treatment of the reprocessing: most 
importantly, constant F was used in computing rhot, while 
the continuum slope evolves during a fiare. Furthermore, 
in our basic model no time delay is assumed in computing 
the reprocessed component, either due to finite light travel 
time or adjustment of hydrostatic equilibrium (Nayakshin 
& Kazanas 2001). However, a possible influen ce of the latter 
effect on the results will be presented in Sec. 4.3.3 in a very 
simple approximation. 

Finally, the /-resolved spectra are computed according 
to the prescription given by Revnivtsev et al. (1999). The 
normalized power spectral density (PSD) is computed as. 



J 



0, 



,iV-l 



JV-1 

E 

fc=0 



Cfc e 



(13) 



(14) 



where the discrete frequencies fj = j/T for j — 
-N/2, . .., N/2; T is the total light curve time, C is the 
mean count rate, while are the number of counts in fc-th 
time bin. The /-resolved spectrum at energy Ei and Fourier 
frequency fj is then defined as 



T 



-A/, 



(15) 



All computations are done for the central black hole 
mass M = 10 M0, accretion rate M — O.OlMEdd, where 
AfEdd ~ L-Edd/{£c^) and e — 0.057 is the accretion efficiency. 
A fraction ^ = 0.2 of total gravitational energy is assumed 
to be dissipated in the disc, with the rest converted to hard 
X-rays (e.g. Di Salvo et al. 2001). 



3 TEST RUNS 

Before presenting the results of the full model computations 
let us consider simpler test runs, which turn out to be helpful 
in understanding the interpretation of the Fourier-/ resolved 
spectra. 

The first, trivial test run is to assume that (1) there is no 
spectral evolution during a flare, i.e. the feedback between 
the soft and hard X-rays is switched off, and (2) all flares 
have the same spectra (no radial dependence of any of the 
parameters). The reprocessed component is not computed. 
All frequency-resolved spectra are then expected to be the 
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Energy (keV) Energy (keV) Energy (keV) 

Figure 2. Test cases of Fourier-/ resolved spectra, plotted as ratio to E~^''^^ power law. Panel (a) sliows an artificial case of a power 
law dependence of r(r) such that r(rnjin) = 1-5 while r(rmax) = 2, i.e. the shorter fiares have harder spectra. The dotted lines show the 
limiting functions E~^'^ and . Labels denote the Fourier frequency with the / < 3 denoting: 1, 0.6, 0.3 and 0.03 Hz. The frequencies 
in all panels are the same. Panels (b) and (c) show the model of PF99 (i.e. no radial dependence of flares' parameters, with spectral 
evolution during fiares), for the PF99 and double exponential heating functions, respectively (see eq. ^ and|2|. Spectra were rescaled 
for clarity. Note opposite behaviour of the spectra for the two li-^(t), and no dependence on / in the range / > 3 Hz, where the PSD is 
formed by superposition of individual flares. 



same and equal to the assumed spectrum, as is indeed the 
case in our computations. 

Second, we set up a run with an artificial dependence 
of spectrum on the flare time scale, r, but still no spectral 
evolution during flares. Specifically, we assume a power law 
r(r) function such, that the shortest fiares (i.e. located at 
r = rinin) have hard spectra, V = 1.5, while the longest ones 
(at r = Tmax) have F = 2. The reprocessed component is not 
computed. Result plotted in Figure ^ show, as expected, 
the high frequency spectra are hard and the spectra soften 
with decreasing /. Below / ~ 3 Hz, where PSD is formed 
by fiares' avalanches, the slope of the energy spectra is no 
longer dependent on /. 

Third, we introduce spectral evolution during flares, but 
neither the radial dependences nor reprocessing. The model 
is thus equivalent to that of PF99, with the same values of 
all parameters, A = 45; /i = 0.8; a = 7; p = 1; Tmin = 1 ms; 
Tmax = 0.15 s, Hq — 0.5, Df) — 0.5. The frequency-resolved 
spectra now too show some Fourier frequency dependence, 
even though all flares have the same time evolution (in time 
units of r). In Fig. |2|b,c the spectra are plotted for the two 
flare profiles (eq. |l| and eq. |2|). The dependencies on / are 
opposite for the two profiles: for the PF99 profile the spectra 
are harder at lower frequencies, while the opposite is true 
for the double exponential profile. All /-resolved spectra are 
identical in the frequency range where individual flares form 
the PSD, i.e. / > 3 Hz (see flg. 2b in PF99), the latter 
value being related to Tmax and chosen to match the break 
in observed PSD. 

The opposite behaviour of energy spectra for / < 3 Hz 
is related to differences in spectral evolution during a flare 
for the two proflles. For the PF99 profile the rise time at soft 
X-rays is shorter than at hard X-rays (see fig. lb in PF99), 
therefore the soft X-ray PSD breaks at higher frequency 
than the hard X-ray PSD. The high-/ energy spectra are 



thus softer. The opposite is true for the double exponential 
profile: the soft X-ray rise time is longer here, the soft X-ray 
PSD breaks at lower /, and the higher-/ energy spectra are 
harder. 



4 FULL MODEL RUNS 

The main qualitative results of the computations with 
the hot skin can actually be easily predicted. Assuming 
Rfi!iic(T) = const and Zh,max('') = const we obtain that the 
maximum illuminating flux is independent of radius. Since 
the gravity decreases with radius, the thickness of the hot 
skin will increase with radius. This is opposite to what is 
observed in usual computations (Nayakshin 2000; Zycki & 
Rozanska 2001), where the illuminating flux is proportional 
to the locally dissipated gravitational energy. Moreover, the 
total luminosity of a flare is proportional to its time scale, 
i'flarc ~ /q i'flarc (i) dt (X T foT both flare profiles. Thus the 
emitted radiation is dominated by emission from outer radii 
(recall that by assumption r oc r^^^) of the emitting region, 
i.e. where the hot skin is thickest. 

4.1 Basic model 

In our basic model we assume the double exponential pro- 
file of the flare heating function, eq. ^, with the following 
parameters: A = 60, fi — 0.7, q = 5, p = 1, Tmin = 1 ms, 
Tmax = 0.15 s. The value of Ho (cxvertical velocity of the 
flare) is such that the height of the flare above the disc is 
equal to its radius at the peak of energy dissipation, while 
Do — 0.7, so that the average spectrum has slope F ~ 1.7. 

The full numerical simulations confirm the above quali- 
tative estimates. The fiux-weighted average thickness of the 
hot skin is plotted in Figure bla, together with the average 
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Figure 3. Results of computations for the basic model witli the double exponential (h(*) (^l- y)- Panel (a) shows the flux-average Tjiot 
and r as functions of radius. The hot skin is thickest at the maximum radius and, as a consequence, the spectrum is hardest there. This 
means that the high-/ spectra, originating at small r, are soft, reversing the trend of hardening with /, as shown in panel (b) — compare 
with Fig. The spectra plotted in panel (b) are for / = 0.03, 0.6, 1, 3, 22 and 40 Hz, ordered by increasing normalization at 2 keV. 
Panel (c) plots the R-F correlation with the Fourier frequency as the control parameter. The correlation holds only for high frequency, 
/ > 3 Hz (solid squares), i.e. when the PSD is formed by individual flares, but it breaks down for avalanches of flares (/ < 3 Hz, empty 
circles). The EW of the Fe Ka line (panel d) shows a complicated behaviour: it follows the R-F correlation for / > 3 Hz; below this 
frequency it is roughly constant, with a weak minimum / ~ 1 Hz. 



slope of the Comptonized continuum. The thickness reaches 
maximum at the maximum radius (recall that the existence 
of the maximum radius comes from the assumption of a 
unique relation between the flare time scale and radius). 
The Comptonized continuum originating at larger distances 
is then indeed harder. Spectral index, F, is smaller than in 
the test run (for the same values of other parameters) be- 
cause of the reduction of the soft reprocessed flux due to the 
hot skin. 

The /-resolved spectra are presented in Figure The 
high frequency, / > /break (/break ~ 2 Hz), spectra are now 
split due to the radial dependence of Thot (compare with 
Fig. ^). Above /break, the higher the frequency the softer 
(intrinsically) the spectra and the larger the amplitude of the 
reprocessed component (so the hardening above ~ 10 keV is 
more significant). The R-F correlation holds for these high 
Fourier frequencies, as shown in Figure The values of R 
and r plotted there at each / were obtained by construct- 
ing an e.g. 'T-spectrum" , ^(r,*) = /h(t) x G[F - Fo(t)] at 



each point in time, where G{x) is a (narrow) gaussian func- 
tion and Fo is the current value of the spectral index. The 
F-spectrum thus describes the current distribution of the 
spectral index, analogously to the energy spectrum describ- 
ing the distribution of energy. Obviously, since F has a well 
defined value at each time step, its distribution is strictly 
speaking a Dirac (5-function. Thus obtained sequence of spec- 
tra can be analyzed analogously to the usual energy spectra. 
In particular /-resolved F-spectra were found, and the mean 
values of F and R from these latter spectra are plotted in 
Fig. An alternative procedure of determining F and R for 
each / would be to fit a reflection model to a spectrum for 
a particular Fourier component. The result would however 
be dependent on the reflection model used, hence the above, 
model independent procedure was adopted. Variation of F 
below / ~ 3 Hz is due to spectral evolution during a flare 
(Fig. B) rather than the presence of the hot skin and, as 
a consequence, the R-T correlation does not hold for those 
lower frequencies. This can also be seen from Fig. Hd, where 
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Figure 4. PSD from the 3 keV lightcurve in tlie model witii radial 
dependence of max- 'h max(^) oc r~^'^. The diameter of the 
flare is assumed constant. The radial dependence of the peak flare 
luminosity influences the distribution of number of events over 
luminosity and, consequently, affects the PSD. In this example 
the PSD has no high-/ break, which is commonly observed in 
black hole binary systems at / ~ 1 Hz. 



the EW of the Fe line is plotted as a function of /. There 
is a clear monotonic dependence EW(/) for / > /break, but 
not so below this frequency. 



4.2 Comparison with observations 

The trends obtained in the presence of the hot skin are 
clearly opposite to what was found from observations of 
black hole binaries. For both Cyg X-1 and GX 339-4 the 
/-resolved spectra are harder when / increases (Revnivtsev 
et al. 1999, 2001). This holds best for the frequency range 
« 0.5-12 Hz, where both single flares and flares' avalanches 
form the PSD. The r(/) dependence seems to flatten both at 
low, / < 0.5 Hz, and high, / > 10 Hz, frequencies (Revnivt- 
sev et al. 1999). Similar conclusion, that the shots contribut- 
ing to low-/ part of PSD have softer spectra than average, 
was reached by Negoro, Kitamoto & Mineshige (2001), us- 
ing the technique of superposed shots proflle. The EW of 
the Fe Kc line decreases monotonically with / for / > 0.5 
Hz, being ^constant below that frequency (Revnivtsev et 
al. 1999, 2001). Again this is incompatible with the model, 
where a rise of EW with /, for / > 1 Hz is observed. At 
low frequency, / < 1 Hz, the model EW(/) is ^constant at 
!=s 40 eV, i.e. the reflection is least efBcient. The observed 
value is « 100 eV, compatible with amplitude of reflection 
R^l (Revnivtsev et al. 1999). 



4.3 Variations on the basic model 

4-3.1 Radial dependence of the flare peak luminosity 

From the discussion above it is clear that the opposite be- 
haviour of Comptonized spectra with Fourier frequency to 
that observed is partly due to the peak flares luminosity 
being constant with radius. Additionaly, the increase of 
luminosity with radius is obviously in disagreement with 
the radial dependence of the gravitational energy release. 



dL/dr{r) oc rr~^. It might therefore seem appropriate to 
modify the Zh,max('') dependence assuming it should decrease 
with radius, so that the longer flares have lower peak lumi- 
nosities. However, any such modification will affect the PSD. 
This is because the PSD is related to the distribution of 
the number of events (individual flares and avalanches) over 
their (total) luminosity. The behaviour is a direct result 
of superposing many small events, fewer bigger events, still 
fewer even bigger events, etc. (e.g. Bak 1997). The break 
at / ~ 2 Hz in PSD in the original model of PF99 seems 
to be related to a break in the distribution dN/dL in this 
model which occurs at L — Tmax^h.max X 27?flarefTieC^/crT (i.e. 
total luminosity of longest flares). Decreasing Zh,max at the 
maximum radius (i.e. longest flares) decreases the range of 
L spanned by individual flares. The resulting distribution 
has no clear break and, as a result, the PSD has no break 
in the whole range of frequencies of interest (Fig. ^). This 
is obviously incompatible with the observed PSD of X-ray 
lightcurves, which universally show breaks at a few Hz. The 
conclusion about lack of break in PSD seems to be indepen- 
dent of t he value of the exponent in the P{r) probability 
(see Sec. 2.1), provided that Zh,max('") is adjusted to give an 
assumed dL/dr(r) (e.g. oc r~^) dependence. Since the total 
emission from each flare is oc r^h.max, the radial distribution 
of emission is dL/dr oc dN/dr x rZh, max oc P(T)rZh.max- For 
the assumed P{t) oc and T(r) oc r^^^ (Eq. the dis- 
tribution becomes dL/dr oc 7'^/'^Zh,max(7')- The PSD plotted 
in Fig. ^ was computed assuming /h,max(7') oc r'^'"\ which 
gives dL/dr{r) oc r~^. 

The effect of the hot skin in this model is insignificant. 
This is because, adjusting the normalization of Zh,max('") so 
that the mean luminosity-weighted Ih ~ 10, the simple time 
average of Zh is much smaller than before. This requires 
larger sizes of flares, -Rflaro, in order for the total luminos- 
ity to remain constant. Hence, the maximum illuminating 
flux, which scales as F oc I/flaro/Zi^ oc l/R^^^.^ (assuming the 
height of maximum dissipation h « -Rflare), is smaller than 
in the model with lh,ma,x{r) — const. As a consequence, no 
correlations is obtained between R and P. The /-resolved 
spectra in this model show slight hardening with /, due to 
flare spectral evolution (see Fig. ^), but the reprocessed 
component is constant with /. 



4-3.2 Modified radial distribution of the stimulated flares 

The Fourier frequency dependence of the spectra becomes 
less signiflcant, when the assumption is made that the spa- 
tial positions of all stimulate d fl ares coincide with that of 

their stimulating flare (Sec 2.2; the distribution 
-1 



of time 

scales P{r) oc is retained for all flares). This is to be 
expected, since the correspondence of variability time scale 
(i.e. /) and physical conditions (i.e. r) is no longer unique. 
In particular, the high-/ spectra are almost identical in this 
case, while the lower-/ ones are somewhat softer compared 
to results of the basic model (Fig. The EW of the Ka 
line shows no frequency dependence being approximately 
constant at ~ 60 eV. 

4-3-3 Delayed response of the reflecting medium 

The properties of the reprocessed component may not be 
always given by the instanteneous illuminating flux, as the 



© 0000 RAS, MNRAS 000, 000-000 



8 P. T. Zycki 



o 



1.2 



1.1 



o 



Id 



11 1 
- r ,= 1.7 

ret 


3-40 Hz /. y - 


- 


/• _ 

jf /l Hz ■ 


1 


1 

\ J* r — 




\ J- 0.03-0.6 Hz - 


, , 1 





5 

Energy 



10 

;keV) 



Figure 5. Fourier /-resolved spectra in the modified model, as- 
suming tliat the locatio n of s timulated flares is the same as their 
stimulating flare (Sec. 4.3.2), i.e. whole avalanche developes in 



the same place. The high-/ (/ > 3 Hz) spectra are now the same, 
since the correspondence of time scale and physical conditions is 
lost. 



accretion disc needs to re-adjust to the condition of hy- 
drostatic equilibrium with changed X-ray radiation pres- 
sure and the Compton temperature (Nayakshin & Kazanas 
2001). Parameters of the hot skin are thus not given by 
the instantaneous flux, if the flux fluctuates rapidly. The 
thickness of the hot skin re-adjusts on the dynamical time 



scale tr. 



tdyn = ' whilc its temperature on the 



thermal time scale, which is usually much shorter (Nayak- 
shin & Kazanas 2001). This effect may influence our results, 
since the response time is approximately equal the flare time 



scale, 



tdy 



Since proper computations of 



time-dependent response of an accretion disc to changing il- 
lumination are not possible at present, we model here the 
response in a rather approximate way, motivated by the re- 
sults of Nayakshin & Kazanas (2001). We simply assume 
that for each flare rhot(i) is determined by illumination flux 



at an eariler time t-i = t — tr 



t — tdyn, i.e. the complex 



response of reflection is approximated as a simple time shift. 
In particular, at the beginning of a flare, for t < idyn, there 
is no hot skin and the reprocessed component is produced 
with amplitude 1. 

The radial dependence of rhot in this model is quali- 
tatively the same as in the basic model, i.e. Thot increases 
with r. The values of rhot('") are now smaller, though, since 
by assumption Thot is fixed at zero for a certain time « r 
during each flare. The radial dependence of F is somewhat 
weaker but the overall trends are virtually the same as in 
our basic model (Fig. ^). With increasing / the spectra 
become harder flrst, then, above / ~ 3 Hz, they soften. The 
EW of the Kq hne has the same /-dependence except for 
a slight overall increase by 15-20 eV. There results are not 
unexpected since the delay is relatively short compared to 
the duration of a flare. With increasing delay, still the same 
qualitative results are obtained. The influence of the hot skin 
becomes less significant but the radial trends are the same, 
if the delay time scale is related to the fiare time scale. In 



the limit of a very long delay the presence of the hot skin 
becomes unimportant, and the model reduces to the original 
model of PF99. 



5 DISCUSSION 

A number of geometrical/physical models of accretion onto 
compact sources can be constructed, if only a part of 
the available information on X-ray emission is used. Sup- 
plementing the usually employed energy spectral analysis 
with time variability and, in particular, correlated spectral- 
timing behaviour, one can attempt to distinguish between 
the models. In this paper predictions for Fourier-resolved 
spectroscopy are presented for the model of magnetic fiares 
avalanches (PF99) above an accretion disc with hot ionized 
skin. 

The basic model formulated by PF99 is quite successful 
in explaining both energy spectra and a number of timing 
characteristics: PSD, hard X-ray time lags and the coher- 
ence function. With the suitably chosen flare proflle (e.g. a 
double exponential with trise ^ tdocay), the auto- and cross- 
correlation function can also be reproduced (Maccarone et 
al. 2000). The dependence of energy spectra on Fourier fre- 
quency, resulting from the spectral evolution during flares 
alone, is in agreement with observations for e.g. the double 
exponential flare proflle, but not for the proflle assumed in 
PF99 (Sec|). 

Our extension of the PF99 model is to assume that there 
exists at least partial correspondence between location of 
a flare and its time scale. The simplest assumption about 
such correspondence is that it follows from projecting the 
flares time scales onto the (inverse of) Keplerian frequency. 
The shortest flares are then located near the inner radius 
of the disc, the longer a flare the farther away it is located. 
Now, in the basic model both the peak heating rate (com- 
pactness parameter) and flare diameter are independent of 
radius, which means that the peak flux of the disc illuminat- 
ing X-rays is independent of radius, too. Since the vertical 
component of gravity decreases with radius (g^ oc {z/r)r~^), 
the thickness of the hot ionized skin increases with radius. 
The fast flares, located where the hot skin is thinest, have 
then softer spectra than the slower flares located where the 
thick hot skin decreases the effectiveness of thermalization of 
X-rays and, consequently, the flux of soft photons for Comp- 
tonization. As a result, the energy spectra corresponding to 
high Fourier frequencies are softer than the lower / spectra, 
for the frequency where the properties of individual flares 
determine the PSD, / > /break, where /break = 1-3 Hz in 
black hole binaries. This is opposite to the observed trends, 
which show steady hardening of the spectra with increasing 
Fourier frequency, at least up to ~ 12 Hz (Revnivtsev et al. 
1999). 

As a consequence, the properties of the reprocessed 
component in the presented model, in particular the Fe Ka 
line, do not agree with the observed trends, either. At high 
/ > /break the reflection amplitude and EW of the Kq line 
increase sharply with /, as a consequence of the inwards- 
decreasing thickness of the hot skin. At low / < 1 Hz the 
model R and EW of the line are roughly constant, how- 
ever R is significantly less than 1 and, consequently, EW is 
smaller than the corresponding value of ~ 100 eV. This is 
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contrary to the observed i? ~ 1 and EW~ 100 eV at these 
low frequencies. 

These results hold true ev en w hen some of the model 
assumptions are modified (Sec. |4.3| ). None of the modifica- 
tions made: the different radial distribution of the stimulated 
flares, radial dependence of the peak flare luminosity or the 
delayed response of reflection, changed the trends of spectral 
parameters with Fourier frequency. It is worth emphasizing 
that the perhaps most obvious of the above modifications 
- the radial distribution of fiares peak luminosity such that 
the total fiares luminosity decreases with radius - not only 
does not help to reproduce the observed EW(/) and r(/), 
but it also changes the PSD so that it no longer agrees with 
observations. It remains to be seen whether a more complex 
formulation of the model is possible, which would simulta- 
neously predict the correct PSD and thickness of the hot 
skin decreasing with radius. 

The same technique of analysis of /-resolved spectra 



can be lapplied to other models mentioned m Sec. I . In the 
dynamic corona model with bulk outflow of the plasma, the 
plasma velocity, /3, is the control parameter of the R-T cor- 
relation (Beloborodov 1999a, b) . The observed /-dependence 
of the energy spectra would require higher /3 closer to the 
center. This is indeed expected, if the X-ray luminosity 
scales with local gravitational energy dissipation. However, 
if the X-ray emission is determined from the properties of 
the flares, assuming their radial distribution as adopted in 
this paper (Sec. Eq. then stronger X-ray emission 
is expected from larger distances, exactly as in the com- 
putations of the hot skin model. It is then likely that (3 
will increase with radius, analogously to the result that the 
thickness of the hot skin increases with radius in our model. 
The result for the /-resolved spectra will then be the same 
as for the model with hot skin. In the alternative geometry 
(the optically thick flow distrupted at certain radius and re- 
placed with some form of hot optically thin flow) , the model 
of Bottcher & Liang (1999) can be used to compute the 
/-resolved spectra. Alternatively, the idea of Kotov, Chu- 
razov & Gilfanov (2001; see also Nowak et al. 1999), that 
the X-ray emitting perturbations are propagating inwards, 
if developed further into a more physical model, could be 
used to predict the /-resolved spectra. 



6 CONCLUSIONS 

• Fourier-resolved X-ray spectroscopy is a powerful tool 
for discriminating various models of accretion. 

• The model of (avalanches of) magnetic flares, radially 
distributed above an accretion disc with hot ionized skin, 
predicts dependencies of various spectral parameters (F, R, 
EW of the Fe Kq line) on Fourier frequency opposite to 
those observed. 
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